INVERSE PROBLEMS IN SPECTRAL GEOMETRY 

KIRIL DATCHEV AND HAMID HEZARI 



I. Introduction 

In 1966, Marc Kac in his famous paper [Kac66] raised the following question: Let VL C IR 2 
be a bounded domain and let 

< A < Ai < A 2 . . . 

be the eigenvalues of the nonnegative Euclidean Laplacian with either Dirichlet or Neu- 
mann boundary conditions. Is Q determined up to isometries from the sequence Ao, Ai . . . ? 
We can ask the same question about bounded domains in IR n , and below we will discuss 
other generalizations as well. Physically, one motivation for this problem is identifying dis- 
tant physical objects, such as stars or atoms, from the light or sound they emit. These inverse 
spectral problems, as some engineers in [Reu07, RWSN09, PeWoRe07, ReWoPe07] have re- 
cently proposed, may also have some interesting applications in shape-matching, copyright 
and medical shape analysis. 

The only domains in W 1 which are known to be spectrally distinguishable from all other 
domains are balls. It is not even known whether or not ellipses are spectrally rigid, i.e. 
whether or not any continuous family of domains containing an ellipse and having the same 
spectrum as that ellipse is necessarily trivial. We can go further and ask the same question 
about a compact Riemannian manifold (M, g) (with or without boundary); can we determine 
(M,g) up to isometries from the spectrum of the Laplace-Beltrami operator A 9 ? Or in 
general, what can we hear from the spectrum? For example, can we hear the area (volume 
in higher dimensions or in the case of Riemannian manifolds) or the perimeter of the domain? 
For the sake of brevity we only mention the historical background for the case of domains. 
In 1910, Lorentz gave a series of physics lectures in Gottingen, and he conjectured that the 
asymptotics of the counting function of the eigenvalues are given by: 

N(X) = H{A,; A, < A} = + 0(V\). 

This asymptotic in particular implies that Area(f2) is a spectral invariant. Hilbert thought 
this conjecture would not be proven in his lifetime, but less than two years later Hermann 
Weyl proved it using the theory of integral equations taught to him by Hilbert. In 1954, 
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Pleijel [Ple54] proved that one knows the perimeter of Q, and in [Kac66], Kac rephrased 
these results in terms of asymptotics of the heat trace 

oo 
3=0 

where the first coefficient a gives the area and the second coefficient gives the perimeter. In 
1967, McKean and Singer [McSi67] proved Pleijel's conjecture that the Euler characteristic 
is also a spectral invariant (this is in fact given by 02) and hence the number of 
holes is known. In 1991, Gordon, Webb and Wolpert [GoWeWo92], found examples of pairs 
of distinct plane domains with the same spectrum. However, their examples were non- 
convex and non-smooth, and it remains an open question to prove that convex domains are 
determined by the spectrum (although there are higher dimensional counterexamples for this 
by Gordon- Webb [GoWe94]) or that smooth domains are determined by the spectrum. 

In this survey we review positive inverse spectral and inverse resonant results for the fol- 
lowing kinds of problems: Laplacians on bounded domains, Laplace-Beltrami operators on 
compact manifolds, Schrodinger operators, Laplacians on exterior domains, and Laplacians 
on manifolds which are hyperbolic near infinity. We also recommend the recent survey of 
Zelditch [Zel04b]. For negative results (counterexamples) we refer the reader to the surveys 
of Gordon [GorOO] and Gordon- Perry-Schueth [GoPeSc05]. 

In the next two sections of the paper we review uniqueness results for radial problems §2, 
and for real analytic and symmetric problems §3. In the first case the object to be identified 
satisfies very strong assumptions (radialness includes full symmetry as well as analyticity) 
but it is identified in a broad class of objects. In this case the first few heat invariants, 
together with an isoperimetric or isoperimetric-type inequality, often suffice. In the second 
case the assumptions on the object to be identified are somewhat weaker (only analyticity 
and finitely many reflection symmetries are assumed) but the identification is only within 
a class of objects which also satisfies the same assumptions, and generic nondegeneracy 
assumptions are also needed. These proofs are based on wave trace invariants corresponding 
to a single nondegenerate simple periodic orbit and its iterations. 

In §4 we consider rigidity and local uniqueness results, where it is shown in the first case that 
isospectral deformations of a given object are necessarily trivial, and in the second case that 
a given object is determined by its spectrum among objects which are nearby in a suitable 
sense. Here the objects to be determined are more general than in the cases considered in 
§2, but less general than those in §3: they are ellipses, spheres, flat manifolds (which have 
completely integrable dynamics), and manifolds of constant negative curvature (which have 
chaotic dynamics). The proofs use these special features of the classical dynamics. 

In §5 we consider compactness results, where it is shown that certain isospectral families are 
compact in a suitable topology. These proofs are based on heat trace invariants and on the 
determinant of the Laplacian, and much more general assumptions are possible than in the 
previous cases. 
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Finally, in §6 we review the trace invariants used for the positive results in the previous 
sections, and give examples of their limitations, that is to say examples of objects which 
have the same trace invariants but which are not isospectral. At this point we also discuss 
the history of these invariants, going back to the seminal paper of Selberg [Sel56]. 

We end the introduction by presenting the four basic settings we consider in this survey: 

1.1. Dirichlet and Neumann Laplacians on bounded domains in IR n . Let f2 be a 
bounded open set with piecewise smooth boundary. Let A^ be the nonegative Laplacian on 
Q with Dirichlet or Neumann boundary conditions. Let 

spec(A^) = (A,)~ , A < Ai < A 2 < • • • 

be the eigenvalues included according to multiplicity, and Uj the corresponding eigenfunc- 
tions, that is to say 

A^Uj = XjUj. 

Recall that A > in the Dirichlet case and A = in the Neumann case. 

1.2. Laplace-Beltrami operators on compact manifolds. Let (M,g) be a compact 
Riemannian manifold without boundary. Let A g = — div 9 grad 9 be the nonnegative Laplace- 
Beltrami operator on M, which we also call the Laplacian for short. Let 

spec(A 9 ) = (Aj)£o, = A < Ai < A 2 < • • • 

be the eigenvalues included according to multiplicity, and Uj the corresponding eigenfunc- 
tions, that is to say 

AgUj = XjUj. 

1.3. Nonsemiclassical and semiclassical Schrodinger operators on W 1 . Let 

V e C°°(R n ;R), lim V(x) = oo, (1.1) 

\x\— >OD 

and let A be the nonnegative Laplacian on W a . Let 

P Vth = h 2 A + V, h> 0, 

Py = Py,l- 

We call Py the nonsemiclassical Schrodinger operator associated to V, and Py^ the semi- 
classical operator. For any h > 0, the spectrum of Py t h on R n is discrete, and we write it 

as 

spec(P^) = (A,)~ o, A < Ai < A 2 < • • • . 

The eigenvalues Xj depend on h, but we do not include this in the notation. We denote by 
Uj the corresponding eigenfunctions (which also depend on h), so that 

Py :hUj = XjUj. 



4 



KIRIL DATCHEV AND HAMID HEZARI 



1.4. Resonance problems for obstacle and potential scattering. In this subsection we 
discuss problems where the spectrum consists of a half line of essential spectrum, together 
with possibly finitely many eigenvalues. In such settings the spectrum contains limited 
information, but one can often define resonances, which supplement the discrete spectral 
data and contain more information. 



1.4.1. Obstacle scattering in M. n . Let O C W 1 be a bounded open set with smooth boundary, 
let il = M. n \ O, and suppose that f2 is connected. Let A^ be the nonnegative Dirichlet or 
Neumann Laplacian on Q. Then the spectrum of Aq is continuous and equal to [0, oo), and 
so it contains no (further) information about Q. One way to reformulate the inverse spectral 
problem in this case is in terms of resonances, which are defined as follows. Introduce a new 
spectral parameter z = V\, with taken so as to map C \ [0, oo) to the upper half plane. 
As Im2 -> + , z 2 approaches [0, oo) and the resolvent (A n — z 2 )^ 1 has no limit as a map 
L 2 {VL) — > L 2 (Q). However, if we restrict the domain of the resolvent and expand the range 
it is possible not only to take the limit but also to take a meromorphic continuation to a 
larger set. More precisely the resolvent 

(where L 2 comp denotes compactly supported L 2 functions and L 2 oc denotes functions which are 
locally L 2 ) continues meromorphically as an operator- valued function of z from {Imz > 0} 
to C when n is odd and to the Riemann surface of log z when n is even. Resonances are 
defined to be the poles of this continuation of the resolvent. Let res(A^) denote the set of 
resonances, included according to multiplicity. See for example [Mel95, Sj602, Zwoll] for 
more information. 



1.4.2. Potential scattering in W 1 . Let Py t h be as before, but instead of (1.1) assume V G 
Co°(IR n ). Then the continuous spectrum of P Vih is equal to [0, oo), but if V is not everywhere 
nonnegative then Py,h m &y have finitely many negative eigenvalues. In either case, the 
resolvent 

(Pv,h - z 2 )' 1 : L 2 comp ->■ L 2 oc , 

has a meromorphic continuation from {Imz > 0} to C when n is odd and to the Riemann 
surface of log z when n is even, and resonances are defined to be the poles of this continuation. 
Let res(iV,h) denote the set of resonances, included according to multiplicity. Again, see for 
example [Mel95, Sj602, Zwoll] for more information. 

1.4.3. Scattering on asymptotically hyperbolic manifolds. The problem of determining a non- 
compact manifold from the scattering resonances of the associated Laplace-Beltrami is in 
general a much more difficult one, but some progress has been made in the asymptotically 
hyperbolic setting. Meromorphic continuation of the resolvent was established by Mazzeo- 
Melrose [MaMe87], and a wave trace formula in the case of surfaces with exact hyperbolic 
ends was found by Guillope-Zworski [GuZw97], which has led to some compactness results: 
see 85.4. 
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2. The radial case 



In this case one makes a strong assumption (radial symmetry) on the object to be spectrally 
determined (whether it is an open set in W 1 , a compact manifold, or a potential) but makes 
almost no assumption on the class of objects within which it is determined. The methods 
involved use the first few heat invariants, and in many cases the isoperimetric inequality or 
an isoperimetric-type inequality. 



2.1. Bounded domains in R n . The oldest inverse spectral results are for radial problems. 
If Q C 1R™ is a bounded open set with smooth boundary, then the spectrum of the Dirichlet 
(or Neumann) Laplacian on Q agrees with the spectrum on the unit ball if and only if fl is 
a translation of this ball. This can be proved in many ways; one way is to use heat trace 
invariants. These are defined to be the coefficients of the asymptotic expansion of the heat 
trace as t — > + : 

oo oo 
3=0 j=0 

where in both the Dirichlet and the Neumann case a is a universal constant times vol(fi), 
and di is a universal constant times vol(<9Q). The left hand side is clearly determined by the 
spectrum, and so the conclusion follows from the isoperimetric inequality. 



2.2. Compact manifolds. Let (M,g) be a smooth Riemannian manifold of dimension n 
without boundary. If n < 6, then the spectrum of the Laplacian on M agrees with the 
spectrum on S n (equipped with the round metric) if and only if M is isometric to S n . This 
was proved by Tanno in [Tan73, Tan80] using the first four coefficients, 00,01,02,03, of the 
heat trace expansion, which in this case takes the form 

00 00 

e~ tx * = Tr e- tA * ~ r n ' 2 ^ aj t j . 

3=0 3=0 

In higher dimensions the analogous result is not known. In [Zel96], Zelditch proves that if 
the multiplicities rrik of the distinct eigenvalues = Eq < E\ < E2 < ■ ■ ■ of the Laplacian 
on M obey the asymptotic = ak 71 ^ 1 + 0(k n ~ 2 ), for some a > as k — > 00 (this is the 
asymptotic behavior for the multiplicities of the eigenvalues of the sphere), then (M,g) is a 
Zoll manifold, that is to say a manifold on which all geodesies are periodic with the same 
period. 



2.3. Schrodinger operators. In general it is impossible to determine a potential V from 
the spectrum of the nonsemiclassical Schrodinger operator A + V. For example, in [McTr81], 
McKean and Trubowitz find an infinite dimensional family of potentials in C°°(R) which are 
isospectral with the harmonic oscillator V(x) = x 2 . 

However, analogous uniqueness results to those above are proved in [DaHeVell] by Ventura 
and the authors, where it is shown that radial, monotonic potentials in IR n (such as for exam- 
ple the harmonic oscillator) are determined by the spectrum of the associated semiclassical 
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Schrodinger operator among all potentials with discrete spectrum. The approach is based 
in part on that of Colin de Verdiere [Col08p], and Guillemin-Wang [GuWa09p] (see also 
[GuStllp, §10.6]), where a one dimensional version of the result is proved. They show that 
an even function (or a suitable noneven function) is determined by its spectrum within the 
class of functions monotonic away from 0. 

The method of proof is similar to that used to prove spectral uniqueness of balls in M. n 
as discussed in §2.1 above. Namely we use the first two trace invariants, this time of 
the semiclassical trace formula of Helffer- Robert [HeRo83], together with the isoperimetric 
inequality. We show that if V, Vq are as in (1.1), if Vo(x) = -R(|x|) where -R(O) = and 
R'(r) > for r > 0, and if spec(P l / i / l ) = spec(P Vih h) up to order 1 o(h 2 ) for h G {hj}JL with 
hj — > + , then V(x) = V (x — x ) for some x G W l . 

The semiclassical trace formula we use is 
Tr(f(P v , h )) = 

jA^(f fm 2 + V)dxdi+^- [ \VV\ 2 f^(\^ + V)dxd^ + 0(h i )), 
where / G Cg°(R). 

Because the spectrum of Py,h is known up to order o{h 2 ) we obtain from (2.2) the two trace 
invariants 

/ dxd£, [ \VV(x)\ 2 dxd£, (2.3) 

J{\tf+V(x)<\} J m 2 +V{x)<x} 

for each A. It follows in particular that V is nonnegative. By integrating in the £ variable, 
we rewrite these invariants as follows: 

f (A - V) n/2 dx, [ \VV(x)\ 2 (X-V) n/2 dx. (2.4) 

J{V(x)<\} J{V(x)<\} 

Using the coarea formula we rewrite the invariants in (2.4) as 

[ X ( I {X ZY^ /2 ds) ds, Hi |W|(A - V) n / 2 ds] ds. 

JO \J{V=s,VV^0} l Vl/ l / J0 \J{V=s} J 

Using the fact that V = s in the inner integrand, the factor of (A — V) n l 2 = (A — s) n l 2 can 
be taken out of the surface integral, leaving 

f (A- s) n/2 h(s)ds, [ (A- s) n/2 I 2 (s)ds, (2.5) 

Jo Jo 

where 

h(s) = [ T^7\dS, Us) = [ \W\dS. (2.6) 

J{V=s,VV^0} I V V I J{V=s} 



^^The implicit rate of convergence here must be uniform on [0, Aq] for each Aq > 0. 



INVERSE PROBLEMS IN SPECTRAL GEOMETRY 7 

We denote the integrals (2.5) by A 1+n / 2 (7i)(A) and A 1+n / 2 (Ii)(X). These are Abel fractional 
integrals of I\ and I 2 (see for example [Zel98b, §5.2] and [GuStllp, (10.45)]), and they can 
be inverted by applying A 1+n / 2 , using the formula 

! A a o-±-A p = 1 „, A a+f) , (2.7) 



r(«) a T(/3) p T(a + /3) 

and differentiating n + 1 times. From this we conclude that the functions I\ and I 2 in (2.6) 
are spectral invariants for every s > 0. 

Integrating I\ and using the coarea formula again we find that the volumes of the sets 
{V < s} are spectral invariants: 

/ h(s')ds' = [ [ -^—dSds'= [ ldx. (2.8) 

JO JO J{V=s',VV^0} |VV I </{V<s} 

From Cauchy-Schwarz and the fact that I\ and I 2 are spectral invariants we obtain 

( I Ids] < ( j^dS [ \VV\dS = [ ^dS [ R'dS, (2.9) 

\J{V=s} J J{V=s} I V V\ J{v=s} J{R=s} R J{R=s} 

when s is not a critical value of V, i.e. by Sard's theorem for almost every s e (0, Ao). On 
the other hand, using the invariants obtained in (2.8) and the fact that the sets {R < s} are 
balls, by the isoperimetric inequality we find 



IdS < / IdS. (2.10) 

{R=s} J{V=s} 



However, 

2 



/ Ids) = [ ]-dS [ R'dS, 

J{R=s} J J{R=s} R J{R=s} 

because l/R' and R' are constant on {R = s}. Consequently 

/ ldS= [ 

JlR=s\ J\\ 



IdS, 

l{R=s} J{V=s} 

and so {V = s} is a sphere for almost every s, because only spheres extremize the isoperi- 
metric inequality. Moreover, 

2 



/ Ids) = [ j^TTjdS [ \VV\dS, 

J{V=s} J J{V=s} I VI/ 1 J{v=s} 



and so IW)" 1 and |VV| are proportional on the surface {V = s} for almost every s, again 
by Cauchy-Schwarz. Using the equation (2.9) to determine the constant of proportionality, 
we find that 

|VVf = R'iR-'is)) 2 = (R^y(s)- 2 = F(s) 
on {V = s}. In other words 

\\7V\ 2 = F(V), (2.11) 

for all x G V^ 1 (s) for almost all s. However, because F(V) ^ when V ^ 0, it follows by 
continuity that this equation holds for all x G l /_1 ((0, oo)). 
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We solve this equation by restricting it to flowlines of W, with initial conditions taken 
on a fixed level set {V = So}, and conclude that, the level surfaces are not only spheres 
(as followed from (2.10)) but are moreover spheres with a common center. Hence, up to 
a translation, V is radial. Since the volumes (2.8) are spectral invariants, it follows that 
V(x) = R(\x\). 

2.4. Resonance problems. We first mention briefly some results for inverse problems for 
resonances for the nonsemiclassical Schrodinger problem when n — 1. In [ZwoOl] Zworski 
proves that a compactly supported even potential V G L X (IR) is determined from the reso- 
nances of Py among other such potentials, and in [Kor05] Korotyaev shows that a potential 
which is not necessarily even is determined by some additional scattering data. 

Analogous results to those discussed in §2.1 hold in the case of obstacle scattering. Hassell 
and Zworski [HaZw99] show that a ball is determined by its Dirichlet resonances among all 
compact obstacles in R 3 . Christiansen [Chr08] extends this result to multiple balls, to higher 
odd dimensions, and to Neumann resonances. As in the other results discussed above, the 
proofs use two trace invariants and isoperimetric-type inequalities, although the invariants 
and inequalities are different here. There is also a large literature of inverse scattering results 
where data other than the resonances are used. A typical datum here is the scattering phase: 
see for example [Mel95, §4.1]. 

In [DaHellp] we prove the analogue for resonances of the result in the previous section 
for semiclassical Schrodinger operators with discrete spectrum. Let n > 1 be odd, and let 
Vo, V G C^°(IR n ; [0, oo)). Suppose V (x) = R(\x\), and R'(r) vanishes only at r = and 
whenever R(r) = 0, and suppose that res(Pv r 0i / l ) = Tes (Pv,h), up to order 2 o(h 2 ), for h G 
{hj}^ =l for some sequence hj — > 0. Then there exists x G K n such that V(x) = V (x — x ). 

Our proof is, as before, based on recovering and analyzing first two integral invariants of 
the Helffer- Robert semiclassical trace formula ([HeRo83, Proposition 5.3], see also [GuStllp, 
§10.5]): 

Tr(f(P v , h ))-f(P 0A ))= (2.12) 
7^ (7 Ml 2 + V) - fm 2 )dxd£ + £ / |W| 2 / (3) (l£| 2 + V)dxdi + <D(h*)) . 

(27Tft,) n \J R 2n 12 J R 2n J 

To express the left hand side of (2.12) in terms of the resonances of Pv,h, we use Mel- 
rose's Poisson formula ([Mel82]), an extension of the formula of Bardos-Guillot-Ralston 
([BaGuRa82]): 

2Tr(cos(t v ^V)-cos(t v ^V)) = ^ e" 1 ^, 0, (2.13) 

Aeres(P v ,h) 

where equality is in the sense of distributions on 1 \ 0. 



The implicit rate of convergence here must be uniform on the disk of radius Aq for each Aq > 0. 
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From (2.13), it follows that if 

g e C °°(R \ 0) is even, (2.14) 

then 

Tr(g(V-h?A + V) - g{\f^WK)) = — V I e'^g^dt. (2.15) 

AGres(P v ,h) 

Now setting the right hand sides of (2.15) and (2.12) equal and taking h — >■ 0, we find that 
/ /(le| 2 + ^)-/(|e| 2 )^, / |V\/| 2 / (3) (le| 2 + ^)^ (2.16) 

JR 2n JR 2n 

are resonant invariants (i.e. are determined by knowledge of the resonances up to o{h 2 )) 
provided that }'{r 2 ) = g{r) for all r and for some g as in (2.14). Taylor expanding, we write 
the first invariant as 

En / f (k) m 2 n I v( X ) k dx+ I YM^L [\i-trf^\\tf+tv(x))dtdxdz. 

k=l K - JM. n Ji™ Jm. 2n m - Jo 

Replacing / by f x , where f x (r) = /(r/A) (note that g x (r) = f x (r 2 ) satisfies (2.14)) gives 

J2 Xn/2 ' k h f / (fc) (l£| 2 K / V(x) k dx + 0(X n/2 ~ m ~ 1 ) 
Taking A — > oo and m-foowe obtain the invariants 



/ f (k) (\en [ v(x) k dx, 

JK" JK. n 



for every k > 1. 



In [DaHellp, Lemma 2.1] it is shown that there exists g satisfying (2.14) such that if f{r 2 ) = 
g(r), then f Rn /«(|£| 2 K ^ 0, provided k > n. 

This shows that 

/ V{x) k dx= [ V (x) k dx (2.17) 
for every k > n, and a similar analysis of the second invariant of (2.16) proves that 

/ V{xf\VV{x)\ 2 dx = [ V (x) k \Wo(x)\ 2 dx (2.18) 

for every k > n. 

We rewrite the invariant (2.17) using V*dx, the pushforward of Lebesgue measure by V, as 

/ V{xfdx= [ s k (V*dx) s =i k VJix {k \{)). (2.19) 

Since V and Vo are both bounded functions, the pushforward measures are compactly sup- 
ported and hence have entire Fourier transforms, and we conclude that 

n— 1 

V*dx = V *dx + ^2 c kbf ] = V *dx + c 5 . 

k=0 
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For the first equality we used the invariants (2.19), and for the second the fact that V*dx is 
a measure. In other words 

vol({V > A}) = vol({T/ > A}) 

whenever A > 0. Moreover, this shows that V*dx is absolutely continuous on (0, oo), and so 
by Sard's lemma the critical set of V is Lebesgue-null on \^ _1 ((0, oo)). As a result we may 
use the coarea formula 3 to write 

V*dx= I {VV^dSds, on(0,oo) 

J{V=s} 

[ |v\/r i rf5 , = f ivvbr 1 ^ 

J{V=s} J{V Q =s} 

for almost every s > 0. Similarly, rewriting the invariants (2.18) as 

/ V{x) k \VV{x)\ 2 dx = [ s k [ \VV\dSds, 



and to conclude that 



we find that 



I \W\dS= [ \W \dS, s>0. 

J{V=s} J{V =s} 



From this point on the proof proceeds as in the previous section. 

To our knowledge it is not known whether such results hold in even dimensions. The higher 
dimensional results discussed above all rely on the Poisson formula (2.13) which is only valid 
for odd dimensions. A similar formula is also true in the obstacle case [Mel83], although 
slightly more care is needed in the definition of cos(t\/A n ) — cos(t^ A K n) because the two 
operators act on different spaces. When n = 1 a stronger trace formula, valid for all t e M, 
is known: see for example [Zwo96, page 3]. When n is even, because the meromorphic 
continuation of the resolvent is not to C but to the Riemann surface of the logarithm, Poisson 
formula? for resonances are more complicated and contain error terms: see [Sj697, Zwo98]. 
A proof based on Sjostrand's local trace formula [Sj697] would be of particular interest, 
firstly because this formula applies in all dimensions and to a very general class of operators, 
and also because it uses only resonances in a sector (and in certain versions, as in [Bon02], 
resonances in a strip) around the real axis. This would strengthen the known results in odd 
dimensions as well as proving results in even dimensions, as one would only have to assume 
that these resonances agreed and not that all resonances do. 



3. The real analytic and symmetric case 



In this case uniqueness results about nonradial objects are obtained, so the assumptions on 
the object to be determined are weaker. However, the assumptions on the class of objects 
within which it is determined are much stronger - in fact they are the same as the assumptions 
on the object to be determined. The two main assumptions are analyticity and symmetry. 

3 If n = 1 we put / |W|-MS = £ x£y - 1(s) \V'(x)\-\ 
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In each case wave invariants are used which are microlocalized near certain periodic orbits, 
as opposed to the non-microlocal heat invariants of the previous section. 



3.1. Bounded domains in IR n . Here the main tool is the following result of Guillemin- 
Melrose [GuMe79b]. When Q C W 1 is a bounded, open set with smooth boundary, they 
prove that Tr(cos(tA/An)) is a tempered distribution in R with the property 



singsuppTr(cos(ty Aq)) C {0}ULsp(fi), 

where Lsp denotes the length spectrum, that is to say the lengths of periodic billiard orbits 
in Q. Moreover, they show that if T e Lsp(f2) is of simple length 4 and 7^ is nondegenerate 5 , 
then for t sufficiently near T we have 



Trcos(t< 



Re 



v/|det(/-P T 



: (t - t + ioy 1 ( 1 + a i(* - T ) j lo s(* - T + i0 ) 



3=1 



+ S(t), 
(3.1) 

where S is smooth near T. Here T" is the primitive length of 7^, which is the length of 
7t without retracing, and a? is the Maslov index of j T (which can be defined geometrically 
but which appears here as the signature of the Hessian in the stationary phase expansion of 
the wave trace). The coefficients cij are known as wave invariants. Viewing the boundary 
locally as the graph of a function /, they are polynomials in the Taylor coefficients of / at 
the reflection points of 7^. In general there is no explicit formula, but they were computed 
by Zelditch [Zel09] in the special case discussed below. 



Assume now n = 2, with coordinates (x,y). Assume further that 



(1) Q is simply connected, symmetric about the x-axis, and dfl is analytic on {y 7^ 0}. 

(2) There is a nondegenerate vertical bouncing ball orbit 7 of length T such that both 
T and 2T are simple lengths in Lsp(fi). 

(3) The endpoints of 7 are not critical points of the curvature of <9f2. 



We recall that a bouncing ball orbit is a 2-link periodic trajectory of the billiard flow, i.e. 
a reversible periodic billiard trajectory that bounces back and forth along a line segment 
orthogonal to the boundary at both endpoints. Without loss of generality we may assume 
that the bouncing ball orbit in assumption (2) is on the y-axis. 

Zelditch [Zel09] proves that if Q and fl' both satisfy these assumptions, and if spec(A^) = 
spec(Afy) (for either Dirichlet or Neumann boundary conditions), then Q = ft' up to a 
reflection about the y-axis. This improves a previous result of Zelditch [ZelOO] where an 
additional symmetry assumption is needed, which in turn improves a previous result of 
Colin de Verdiere [Col84] where rigidity is proved in the class of analytic domains with two 

4 This means that only one periodic orbit (up to time reversal), 7t, has length T. 

5 This means that 77- is transversal to the boundary and Pt the linearized Poincare map of 7t, which is 
the derivative of the first return map, does not have eigenvalue 1. 
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reflection symmetries. Under the above assumptions, for e > sufficiently small, there exists 
a real analytic function /: (—e,e) — > R such that 

ttn{\x\ <e} = {(x,y): \x\ < e, \y\ < f{x)}. 

To prove the theorem it is enough to show that the Taylor coefficients of / at are determined 
by spec(f2) (up to possibly replacing f(x) by f(—x)). Zelditch does this by writing a formula 
for the coefficients a,j of (3.1), which are determined by spec(f2), applied to 7 and to 7 2 (the 
iteration of 7): 

%(7 r ) =Mr)f^\0) + ^(r)/^ +1 )(0)/( 3 )(0) 

+ [terms containing f^ k \0) only for k < j\ . 

Here Aj(r) and Bj(r) are spectral invariants which are determined by the first term of (3.1). 
One can show that (A,(l), Bj(l)) as a vector is linearly independent from (A,-(2), Bj(2)). 
Hence, by an inductive argument, if f^ 3 \0) 7^ 0, all the coefficients are determined (up to a 
sign ambiguity for /' 3 '(0), which corresponds to reflection about the y-axis). The condition 
/^(0) 7^ is equivalent to assumption (3) above, and in [Zel09, §6.9] Zelditch outlines a 
possible proof in the case where /^(0) = 0. 

In [HeZelOp], Zelditch and the second author prove that bounded analytic domains (let" 
with ± reflection symmetries across all coordinate axes, and with one axis height fixed (and 
also satisfying some generic non-degeneracy conditions) are spectrally determined among 
other such domains. This inverse result gives a higher dimensional analogue of the result 
discussed above from [Zel09], but with n axes of symmetry rather than n — 1. To our 
knowledge, it is the first positive higher dimensional inverse spectral result for Euclidean 
domains which is not restricted to balls. The proof is based as before on (3.1) and on 
formulas for the a,j(Y)i but there are additional algebraic and combinatorial complications 
coming from the fact that Taylor coefficients must be recovered corresponding to all possible 
combinations of partial derivatives. These complications are very similar to those that arise 
for higher dimensional semiclassical Schrodinger operators discussed below in §3.3 



3.2. Compact manifolds. To our knowledge all uniqueness results in this category are 
about surfaces of revolution. In [Ber76] and [Gur95], Berard and Gurarie show that the 
joint spectrum of A g and d/d6 (the generator of rotations) of a smooth surface of revolution 
determines the metric among smooth surfaces of revolution, by reducing the problem to a 
semiclassical Schrodinger operator in one dimension. 

In [BrHe84], Bruning-Heintze show that the spectrum of A 9 alone determines the metric of 
a smooth surface of revolution with an up-down symmetry among such surfaces. They prove 
that the spectrum of A g determines the ^-invariant spectrum (but not necessarily the full 
joint spectrum), allowing them to apply a result of Marchenko [Mar52, Theorem 2.3.2] for 
one-dimensional Schrodinger operators. 

In [Zel98b], Zelditch proves that a convex analytic surface of revolution satisfying a non- 
degeneracy condition and a simplicity condition is determined uniquely by the spectrum 
among all such surfaces. He uses analyticity and convexity to show that the spectrum 
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determines the full joint spectrum of A g and d/d6, reducing the problem to a semiclassical 
Schrodinger operator in one dimension. 

3.3. Schrodinger operators. When n — 1, Marchenko [Mar52] shows that an even poten- 
tial is determined by the spectrum of the associated nonsemiclassical Schrodinger operator 
among all even potentials. More specifically, he proves in [Mar52, Theorem 2.3.2] that a 
Schrodinger operator on [0, oo) is determined by knowledge of both the Dirichlet and the 
Neumann spectrum. The result for even potentials on R follows from the result on [0, oo) 
as follows: Let V E C°°(M) obey hin^^oo V(x) = oo. If Uj is the eigenfunction of — + V 
corresponding to the eigenvalue Aj, then Uj has exactly j zeros and they are all simple (see 
[BeSh91, Chapter 2, Theorem 3.5]). If V is even then every eigenfunction is either odd or 
even and this result shows that the parity of Uj is the same as the parity of j. In particular 
(A^)°^ with A^ = \ 2 j is the spectrum of 

d 2 

— — + V on L 2 (\0, oo)) with Neumann boundary condition, 

dx l 

and (A")°^ with A" = A 2 j+i is the spectrum of 
d 2 

— — - + V on L 2 ([0, oo)) with Dirichlet boundary condition. 
dx 2 

This reduces the problem on R to the result of Marchenko. 

However, noneven potentials may have the same spectrum: indeed McKean-Trubowitz 
[McTr81] construct an infinite-dimensional family of potentials which have the same spec- 
trum as the one-dimensional harmonic oscillator V(x) = x 2 . 

In [GuUr07], Guillemin-Uribe consider potentials V in M. n which are analytic and even in all 
variables, which have a unique global minimum U(0) = 0, which obey liminf^i^oo V(x) > 0, 
and such that the square roots of the eigenvalues of Hess V(0) are linearly independent over Q. 
They show that such potentials are determined by their low lying semiclassical eigenvalues, 
that is to say by spec(Pv,h) H [0,e] for any e > 0. In [Hez09], the second author removes 
the symmetry assumption in the case n = 1 but assumes V"'(0) ^ 0, and for n > 2 replaces 
the symmetry assumption by the assumption that V(x) = }\x\ 1 . . . x^) + x^g(x 2 , . . . , xl). 
In [CoGu08p] and [Col08p], Colin de Verdiere-Guillemin and Colin de Verdiere give another 
proof of this result for the case n — 1, and in [GuUrll], Guillemin-Uribe give another proof 
in the higher-dimensional case. 

The proofs in [GuUr07, CoGu08p, Col08p, GuUrll] are based on quantum Birkhoff normal 
forms. These are a quantum version of the Birkhoff normal forms of classical mechanics. 
In the classical case, one constructs a symplectomorphism which puts a Hamiltonian func- 
tion into a canonical form in a neighborhood of a periodic orbit. In the quantum case, 
one constructs a Fourier integral operator associated to this symplectomorphism which puts 
a pseudodifferential operator which is a quantization of this Hamiltonian into a canonical 
form, micro locally near the periodic orbit. Quantum Birkhoff normal forms were developed 
by Sjostrand [Sj692] for semiclassical Schrodinger operators near a global minimum of the 
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potential. Guillemin [Gui96] and Zelditch [Zel97, Zel98a] put the Laplace-Beltrami opera- 
tor on a compact Riemannian manifold into a quantum Birkhoff normal form. Sjostrand- 
Zworski [SjZw02] and Iantchenko-Sjostrand-Zworski [IaSjZw02] studied general semiclassical 
Schrodinger operators on a manifold at nondegenerate energy levels. 

The proof in [Hez09] (that the Taylor coefficients of the potential at the bottom of the well 
are determined by the low-lying eigenvalues) is based on Schrodinger trace invariants. These 
are coefficients of the expansion 



where \ £ Co°0^) is 1 near and is supported in a sufficiently small neighborhood of 0. The 
coefficients a, in dimension n = 1 have exactly the form (3.2) (and in higher dimensions they 
have the same form as the higher dimensional coefficients of the wave trace on a bounded 
domain) and hence, once this fact is established, the remainder of the uniqueness proof is 
the same for both problems. 

3.4. Resonance problems. The case of an analytic obstacle with two mutually symmetric 
connected components is treated by Zelditch [Zel04a] by using the singularities of the wave 
trace generated by the bouncing ball between the two components. In [Zwo07] , Zworski gives 
a general method for reducing inverse problems for resonances on a noncompact space to 
corresponding inverse problems for spectra on a compact space. 

In [Ian08] Iantchenko considers potentials V in W 1 which are analytic and even in all variables, 
which have a unique global maximum at V(0) = E, which extend holomorphically to a sector 
around the real axis and obey liminf|. c |_ >00 V(x) = in that sector, and such that the square 
roots of the eigenvalues of Hessl / (0) are linearly independent over Q. He uses the quantum 
Birkhoff normal form method of [GuUr07] to recover the Taylor coefficients of the potential 
at the maximum and to show that potentials V in this class are determined by the resonances 
in a small neighborhood of E. 



In this section we consider results which show nonexistence of nontrivial isospectral defor- 
mations. 

4.1. Bounded domains in W 1 . In [MaMe82], Marvizi-Melrose introduce new invariants 
for strictly convex bounded domains flcl 2 based on the length spectrum, associated with 
the boundary. They show that, for m G N fixed, 

sup {£(7) : 7 is a periodic billiard orbit with m rotations and n reflections} 



Tr(e- itP ^ h X (P v , h )) = J>(i)^, h + , 



3=0 



4. Rigidity and local uniqueness results 



■oo 




n — > 00, 



k=i 
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where L denotes the length. Then they introduce the following non- coincidence condition 
on Q, which holds for a dense open family (in the C°° topology) of strictly convex domains: 
suppose there exists e > such that if 7 is a closed orbit with L(dQ) — e < L(j) < L(dfl), 
then 7 consists of one rotation. They show that under this condition, the coefficients Cfc >m 
are spectral invariants, and they use the invariants c^i and to construct a two-parameter 
family of planar domains which are locally spectrally unique (meaning that each domain has 
a neighborhood in the C°° topology within which it is determined by its spectrum). The 
two parameter family consists of domains which are defined by elliptic integrals, and which 
resemble ellipses, but which are not ellipses. 

In [GuMe79a], Guillemin-Melrose consider the Laplacian on an ellipse Q given by x 2 /a + 
y 2 /b = 1, with a > b > 0, and with boundary condition given by 

du/dn = Ku, on <9f2, (4.1) 

where K G C°°(dVL) and is even in both x and y. They show K is determined by spec(An^), 
where is the Laplacian on Q with boundary condition (4.1). 

To explain their method, let us introduce some terminology. For T > the length of a 
periodic orbit, the fixed point set of T, denoted by Y T , is the set of (q : rj) G B*dQ, the 
coball bundle of dQ, such that the billiard orbit corresponding to the initial condition (q, 77) 
is periodic and has length T. For a more general domain there will often be only one periodic 
orbit of length T (up to time reversal), but an ellipse, because of the complete integrability of 
its billiard flow, always has one or several one-parameter families of such orbits. Guillemin- 
Melrose prove that, in the case of the ellipse, for any T which is the length of a periodic 
orbit such that L(dQ) — T > is sufficiently small, Yt has one connected component T (up 
to time reversal). This connected component is necessarily a curve which is invariant under 
the billiard map. Moreover, they show that the asymptotic expansion of 

Tr(cos(tv/A^)) - Tr(cos(tv/A^)) 

in fractional powers of t — T has leading coefficient 

f K 

/ -y== 2 d ^- (4.2) 

Jr a/ 1 - T 

Here /xr is the Leray measure on T. They then show that, under the symmetry assumptions, 
K is determined from a sequence of such integrals for Tj with Tj tending to L(d£l) from 
below. 

In [HeZelOp], Zelditch and the second author prove that an ellipse is infinitesimally spectrally 
rigid among C°° domains with the symmetries of the ellipse. This means that if Q Q is an 
ellipse, and if p e is a smooth one-parameter family of smooth functions on dfl which are 
even in x and y, and if Q e is a domain whose boundary is defined by 

dfl € = {z + p e (z)n z : z G dfl }, 

and if spec(f2o) — spec(f2 e ) for e G [0, eo), then the Taylor expansion of p t vanishes at e = 0. 
In particular, if p depends on e analytically, the deformation is constant. The proof uses 
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Hadamard's variational formula for the wave trace: 



d_ 
de 



Trtcos^v^)) 



t 



e=0 



d ni d n2 Sci (t, z, z) 



an 



de 



Pe{z) dz, 



e=0 



(4.3) 



where d ni and <9„ 2 denote normal derivatives in the first and second variables respectively, 
Sq is the kernel of sm(tW Aq )/ ^Aq . They then use (4.3) to prove that for any T in the 
length spectrum of fl , the leading order singularity of the wave trace variation is, 



d 
de 



Tr(cos(tv^)) 



e=0 



FCY T 1 v 



p< 



e=0 



v/l-H 2 d// r )(t-T + iO)-i}, 



(4.4) 



modulo lower order singularities, where the sum is over the connected components T of the 
set Yt of periodic points of the billiard map on B*8Qq (and its powers) of length T, and 
where dfir is as in (4.2). As before, if L(dil ) — T > is sufficiently small, there is only one 
connected component and the sum has only one term. For an isospectral deformation, the 
left hand side of (4.4) vanishes, and hence the integrals 



/ r (C/-) v/r ^* r 



vanish when L(<9fi ) — T > sufficiently small. From this point on proceeding as in 
[GuMe79a] above one can show that | e _ Pe = 0, and reparametrizing the variation one can 
show that all Taylor coefficients of the variation are 0. In [HeZelOp] it is shown that expan- 
sions of the form (4.3) and (4.4) hold more generally and in higher dimensions; indeed (4.3) 
holds for any C 1 variation of any bounded domain, and a version of (4.4) holds whenever 
the fixed point sets Y T are clean. These formulas may be useful for example in a possible 
proof of spectral rigidity of ellipsoids. 



4.2. Compact manifolds. In [Tan80] Tanno, uses heat trace invariants to show local spec- 
tral uniqueness of spheres in all dimensions. This means that there is a C°° neighborhood 
of the round metric on the sphere within which this metric is spectrally determined. In 
[Kuw80], Kuwabara does this for compact flat manifolds, and in [Sha09], Sharafutdinov does 
this for compact manifolds of constant negative curvature. 

In [GuKa80a, GuKa80b], Guillemin-Kazhdan prove that a negatively curved compact man- 
ifold (M, g) with simple length spectrum is spectrally rigid if its sectional curvatures satisfy 
the pinching condition that for every x G M there is A(x) > such that \K/A + 1| < 1/n, 
where K is any sectional curvature at x (note that the pinching condition is satisfied for all 
negatively curved surfaces because in that case there is only one sectional curvature at each 
point x and we may take A(x) = —K). Spectrally rigid here means if g € is a smooth family of 
metrics on M with g = g and with spec(A s J = spec(A s ), then (M, g e ) is isometric to (M, g) 
for every e. They further use a similar method of proof to establish a spectral uniqueness 
result for Schrodinger operators on these manifolds. The pinching condition was relaxed by 
Min-Oo in [Min86] and removed by Croke-Sharafutdinov in [CrSh98], and the result was 
extended to Anosov surfaces with no focal points by Sharafutdinov-Uhlmann in [ShUhOO]. 
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4.3. Schrodinger operators. In [Hezllp], the second author considers anisotropic har- 
monic oscillators: V(x) + • • • + a^Xni where the CI j £1X6 linearly independent over Q. It 
is shown that if V t (x) is a smooth deformation of V(x) within the class of C°° functions which 
are even in each Xj, and if spec(Py ) i) = spec(Py ei i) for all e G [0, e), then the deformation is 
flat at e = 0, just as in the infinitesimal rigidity result for the ellipse in §4.1. 

5. Compactness results 

5.1. Bounded domains in W 1 . In [Mel83p], Melrose uses heat trace invariants and Sobolev 
embedding to prove compactness of isospectral sets of domains O C I 2 in the sense of the 
C°° topology on the curvature functions in C°°(dQ). This result allows the possibilty of a 
sequence of isospectral domains whose curvatures converge but which 'pinch off' in such a 
way that the limit object is not a domain, but in [Mel96] points out that this possibility can 
be ruled out using the fact that the singularity of the wave trace at t — is isolated. 

In [OsPhSa89], Osgood-Phillips-Sarnak give another approach to this problem based on the 
determinant of the Laplacian. This is defined via the analytic continuation of the zeta 
function 

oo 

Z(*) = j>J', detA n = e -*<°>. 
i=i 

They consider the domain Q as the image of the unit disk D under a conformal map F, 
with e^go the induced metric on D, where = log \F'\ is a harmonic function. Thus 
is determined by its boundary values, and the topology of [OsPhSa89] is the C°° topology 
on 4>\dDi an d in this case pinching degenerations are ruled out automatically. In [HaZe99], 
Hassell-Zelditch give a nice review of these results and an application of these methods to 
the compactness problem for isophasal obstacles in IR 2 . 

To our knowledge there is no compactness result in higher dimensions. 

5.2. Compact manifolds. In [OsPhSa88a, OsPhSa88b], Osgood-Phillips-Sarnak extend 
their determinant methods to the case of surfaces and prove that the set of isospectral 
metrics on a given Riemannian surface is sequentially compact in the C°° topology, up to 
isometry. In [ChYa89] and in [BrPeYa89], Chang- Yang and Brooks- Perry- Yang give com- 
pactness results for isospectral metrics in a given conformal class on a three dimensional 
manifold. In [BrPePe92], Brooks-Perry-Petersen prove compactness for isospectral families 
of Riemannian manifolds provided that either the sectional curvatures are all negative, or 
that there is a uniform lower bound on the Ricci curvatures. In [Zho97], Zhou shows that 
on a given manifold, the family of isospectral Riemmanian metrics with uniformly bounded 
curvature is compact, with no restriction on the dimension. 

5.3. Schrodinger operators. In [Brii84], Bruning considers Schrodinger operators A g + V 
on a compact Riemannian manifold (M,g), where V G C°°(M), and proves that if the 
dimension n < 3, then any set of isospectral potentials is compact. In higher dimensions 
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he proves the same result under the additional condition that the H s norm of V for some 
s > 3(n/2) — 2 is known to be bounded by some constant C. In [Don05], Donnelly improves 
this condition to s > (n/2) — 2, and derives alternative compactness criteria: he shows 
that isospectral families of nonnegative potentials are compact in dimensions n < 9. If one 
considers instead Ag + '-fV, he shows that a family of potentials which is isospectral for more 
than (n/2) — 1 different values of 7 is compact. In particular, this implies compactness of 
families which are isospectral for the semiclassical problem h 2 A g + V. 

5.4. Resonance problems. Let (X ,g ) be a conformally compact surface, which is hy- 
perbolic (has constant curvature) outside a given compact set K C X . This means that, if 
K Q is taken sufficiently large, then X \ K is a finite disjoint union of funnel ends, which is 
to say ends of the form 

(0, oo) r x Si, dr 2 + f cosh 2 (r)d6 2 , (5.1) 

where £ 7^ may vary between the funnels. Then the continuous spectrum of A 9o is given by 
[1/4, 00), and the point spectrum is either empty or finite and contained in (0, 1/4) (and there 
is no other spectrum). If we introduce the spectral parameter z = a/A — 1/4, where is 
taken to map C\[0, 00) to the upper half plane, then the resolvent (A 9 — 1/4 — z 2 )^ 1 continues 
meromorphically from {Imz > 0} to C as an operator L^ omp — > L 2 oc . This meromorphic 
continuation can be proved by writing a parametrix in terms of the resolvent of the Laplacian 
on the ends (5.1), which in this case can be written explicitly in terms of special functions: 
see [MaMe87] for the general construction, and [GuZw95, §5] for a simpler version in this 
case. 

Borthwick- Perry [BoPell] use a Poisson formula for resonances due to Guillope-Zworski 
[GuZw97] and a heat trace expansion to show that the set of surfaces which are isoresonant 
with (Xo, go) and for which there is a compact set K C X such that (Xo\Kq, go) is isometric 
to (X \ K, g) is compact in the C°° topology, improving a previous result of Borthwick- 
Judge-Perry [BoJuPe03]. They also prove related but weaker results in higher dimensions. 

6. Trace invariants and their limitations 

6.1. Bounded domains in IR n . For with f2 C K" a bounded smooth domain we have 
seen two kinds of trace invariants. The first are heat trace invariants, which are the coeffi- 
cients cij of the expansion 

00 

Tre- tA « ^r n ' 2 Y,^ 3l \ ^0+, 
3=0 

are given by integrals along the boundary of polynomials in the curvature and its derivatives. 
These are equivalent to the invariants obtained from coefficients of the expansion of the wave 
trace Tr cos(ty/A n ) at t = 0. 

The other kind are wave trace invariants obtained from coefficients of the expansion of the 
wave trace at the length of a periodic billiard orbit, always assumed to be nondegenerate 
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and usually assumed to be simple. In this case the formula, as already mentioned in (3.1), is 



Re 



y/det(I -P T ) \ j^i 



+ S(t), 
(6.i; 



where 77- is the simple periodic orbit of length T, and where the coefficients bj are polynomials 
in the Taylor coefficients at the reflection points of 771 of the function of which the boundary 
is a graph. Because of this requirement on the periodic orbit, positive inverse results of the 
kind described above, which are based on the wave trace, always require generic assumptions 
such as nondegeneracy and simple length spectrum. Although there has been some work 
on the degenerate case, such as [Pop98] by Popov, it does not seem to have led yet to 
uniqueness, rigidity, or compactness results. However, in [MaMe82], Marvizi- Melrose obtain 
information from invariants at lengths approaching the length of dQ (see §4.1 above for more 
information) . 

Another limitation comes from the fact that domains can have the same trace invari- 
ants without being isospectral. That is to say, we can construct Q and Q' such that 
Tr(cos(t v / An)) - Tr(cos(t v / A^)) G C°°(R) (recall that the wave trace invariants are the 
coefficients in the expansion of the wave trace near a singularity, as in (6.1)), but spec(A^) 7^ 
spec(An')- This was done by Fulling- Kuchment in [FuKu05], following a conjecture of 
Zelditch [Zel04b], where the following types of domains are considered (these were first 
introduced by Penrose to study the illumination problem, and then shown by Lifshits to be 
examples of nonisometric domains with the same length spectrum): 





Figure 1. Two domains Q and Q' with Tr(cos(tA/A^)) - Tr(cos(t\/A^7)) G 
C°°(R), but spec(A n ) ^ spec(A n ,). 
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These two domains are obtained by taking a semi-ellipse and adding two asymmetric bumps 
A, B and A', B', with A = A' and B = B', such that the foci are left unperturbed (as in 
the figure). Then one adds bumps C and C, the small bumps in the middle which are in 
between the foci, such that C ^ C but C and C are reflections of one another. These two 
domains are not isometric but have the same heat invariants, because heat invariants are 
given by integrals along the boundary of polynomials in the curvature and its derivatives - 
indeed we have freedom to 'slide' C back and forth along the boundary without changing 
any heat invariants, although this is not the case for the wave trace invariants. 

We now show that Tr(cos(tv^)) - Tr(cos(t v / A^)) e C°°(K). This is because of the 
following separation of the phase spaces 6 B*dQ and B*dQ' into two disconnected rooms 
each, which are invariant under the billiard maps of the domains, and which we denote Ri, 
R 2 , R[, and R' 2 , and which have the property that R\ U R2 — B*dVt and R[ U R' 2 = B*dQ! . 
These are defined as follows: R\ is the set of points in B*dQ whose billiard flowout intersects 
the part of the boundary strictly in between the two foci, R2 is the set of points in B*dQ 
whose billiard flowout intersects the part of the boundary which is strictly outside the two 
foci but on the axis of the ellipse or below (and similarly for R[ and Rl 2 ). These two sets are 
disjoint because billiards in an ellipse which intersect the major axis in between the two foci 
once do so always. Now we make the generic assumption that no trajectory which passes 
through the two foci in the initial semi-ellipse is periodic. Because R\ is isometric to R[, and 
R 2 is isometric to R' 2 , we have Tr(cos(t>/Afi)) — Tr(cos(t^/Sfy)) £ C°°(R). This is because 
the singularities of Tr(cos(tV Ac>)) occur at t — T, where T is the length of a periodic orbit, 
and only depend on the structure of B*dQ in an arbitrarily small neighborhood of the orbits 
of length T. 

To show that spec(A^) 7^ spec(Afy), Fulling-Kuchment use a perturbation argument based 
on Hadamard's variational formula for the ground state to show that, for suitably chosen 
small C, the ground states are not the same. 

6.2. Compact manifolds. As we have already mentioned, heat trace invariants can be 
defined for compact manifolds (M, g) as well. In the boundaryless case the expansion takes 
the form 



where the dj are given by integrals on M of polynomials in the curvature and its derivatives. 
Half powers of t appear only when there is a boundary, as in the case of domains considered 
above. Once again, these invariants are equivalent to the invariants obtained from coefficients 
of the expansion of the wave trace Trcos(tAyA^) at t — 0. In analogy with the example given 
in the previous section, we can construct manifolds which are not isometric but which have 
the same heat invariants by taking a sphere, adding two disjoint bumps, and moving them 
around. For suitable choices of bumps it should be possible to make the length spectra 
nonequal, as a result of which the manifolds will be nonisospectral. It seems to be an 
open problem, however, to find an example of two manifolds (M,g) and (M',g') which are 



'Recall that B*dfl is the ball bundle, the fibers of which are intervals [—1,1]. 
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nonisospectral but which have Tr(cos(t^ lt / r \ j )) — Tr(cos(iy / A^)) G C°°(R), that is to say 
which have identical wave trace invariants. 

In this setting the wave trace expansion was established by Duistermaat-Guillemin [DuGu75] , 
building off of previous work by Colin de Verdiere [Col73] and Chazarain [Cha74]. It is a 
generalization of Selberg's Poisson formula [Sel56] to an arbitrary compact boundaryless 
Riemannian manifold. For T the length of a simple nondegenerate periodic geodesic 7t, it 
takes the form 

rpt / OO \ 

TreV A s =i°r - g-T + iO)- 1 1 + V bAt - TV log(t - T + iO) + S(t), 

V|det(/-P T )| l \ U J 

where S(t) is smooth near T. Using quantum Birkhoff normal forms, in [Zel98a], Zelditch 
shows these coefficients bj to be integrals of polynomials in the metric and its derivatives 
along 7j\ See also [Zel99] for a more detailed survey on wave invariants. Because of this 
very local nature of these invariants, to prove uniqueness results one must either assume 
analyticity (as is done in the results discussed above) or find a way to combine information 
from many different orbits (no one seems to have been able to do this so far). 



6.3. Schrodinger operators. In the setting of semiclassical Schrodinger operators the ana- 
logue of the Duistermaat-Guillemin wave trace is the Gutzwiller trace formula near the 
length T of a periodic trajectory of the Hamiltonian vector field H p in p~ 1 (E), where 
p(x,0 = \Z\ 2 + V(x): 

iS jh 00 

Tre-^ Pv ' h - E ^ h x(P Vh ) ~ Va 77 fr J , a 07 = 5 (t - T), 

for t near T, where \ ^ Co°(^) has \ — \ near E. Here the sum in 7 is over periodic 
trajectories in p _1 (0) of length T and the distributions whose singular support is 

contained in {T}. This formula goes back to work of Gutzwiller [Gut71], and was proved 
in various degrees of generality and with various methods by Guillemin-Uribe [GuUr89] 
and Combescure-Ralston-Robert [CoRaRo99] (see also this last paper for further history 
and references). This formula is also valid for more general pseudodifferential operators 
of real principal type, so long as E is a regular value of the principal symbol p and so 
long as the periodic trajectories in p~ 1 (E) are nondegenerate (so that the determinants in 
the denominator are nonzero). In particular it also applies on manifolds. In [IaSjZw02], 
Iantchenko-Sjostrand-Zworski use it to recover quantum and classical Birkhoff normal forms 
of semiclassical classical Schrodinger operators, at nondegenerate periodic orbits. When the 
energy level is degenerate, the Gutzwiller trace formula becomes more complicated: see §3.3 
for a discussion of the case where p has a unique global minimum at E, and see for example 
[BrPaUr95] and [Khu97] for other cases. 

In [Col08p], Colin de Verdiere gives an example of a pair of potentials V ^ V G C°°(R) 
such that spec(Py t h) — spec{Pv,h) up to 0(h°°), so that in particular all semiclassical trace 
invariants for these two potentials agree. He conjectures however, that the spectra are not 
equal. In [GuHellp], Guillemin and the second author construct a pair of potentials, which 



22 



KIRIL DATCHEV AND HAMID HEZARI 



are perturbations of the harmonic oscillator analogous to the perturbations of the semi- 
ellipse discussed in §6.1, and which have different ground states and hence different spectra, 
although the spectra still agree up to 0(h°°). 
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